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The D234 action for light quarks 
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We investigate a new light fermion action (the "D234" action), which is accurate up to ©(a^) and tadpole- 
improved 0{aas) errors. Using D234 with Symanzik- and tadpole-improved glue we find evidence that continuum 
results for the quenched hadron spectrum (pion, rho and nucleon) can be obtained on coarse lattices. 



1. INTRODUCTION 

It has recently become clear that improved ac- 
tions [ |l|, 1^ are a powerful and practical tool 
for obtaining continuum results from the lattice. 
This has been demonstrated for the T with the 
non-relativistic QCD action [ ^, and also for the 
static quark potential and charmonium with im- 
proved glue [ ^. In the latter case, continuum 
results were obtained from lattices as coarse as 
a = 0.4 fm. 

At tree level, improved glue has 0{a'^) er- 
rors, much smaller than the 0{a^) errors of 
the Sheikholeslami-Wohlert (SW) action [ D for 
fermions, which is the most improved light quark 
action used in simulations up to now (see other 
contributions at this conference). In order to 
bring the quark action up to the same level of ac- 
curacy as the glue, we have started a program of 
investigating improved light fermion actions. We 
here report results for the "D234" action in the 
quenched approximation. The D234 action uses 
a cubic derivative term to improve the naive lat- 
tice Dirac operator. It also contains second and 
fourth order derivative terms (as well as the a ■ F 
"clover" term) with fine tuned coefficients to elim- 
inate the "spatial doublers" at large momentum 
and minimize the number of "time doublers" at 
low momentum without introducing errors larger 
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than 0{a^). 

As with other improved actions, we find that 
tadpole improvement [ |^ of the coefficients is 
essential in obtaining accurate results with the 
D234 action. At tree level it can be implemented 
by simply replacing each gauge link U in the ac- 
tion hy U = U /uq, where (wq)^ is the expectation 
value of the plaquette. 

2. D234 ACTION 

The D234 fermionic operator used in the sim- 
ulations described here is: 



mo -I 
ra 



5: 7.AW-6a27.AWA(2) 



(1) 



where h = 1/6, r = 2/3, c= 1/12 



a 

+ U-f,{x)'iJj{x ~ - 2il;{x) 



(2) 



(3) 



and for we use the clover representation [ ^ . 

Classically the D234 action has 0{a^) errors. 
One can show this by using a change of variable 
to relate it to 7pA|i^'(l - ^A^^^), which differs 
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at 0{a'^) from the continuum Dirac operator. For 
the rest of this paper we ignore the change of vari- 
able since it does not affect spectral quantities. 
At the 1-loop level there are additional 0{ag^) 
errors. The coefficients r and c are chosen so that 
there is only one low momentum doubler in the 
free dispersion relation, as can be seen in fig. |l|. 
For the glue we used the Symanzik- and tadpole- 
improved action described in [ which has 
0{a^) and 0{a'^g^) errors. 




Figure 1. Massless dispersion relations for free 
D234 (sohd hnes), Wilson or SW (dotted hne) 
and continuum (dashed line) fermions. Momen- 
tum p is along the (1,1,0) direction. Beyond the 
D234 branch point, the real part of the two con- 
jugate roots is shown. 



3. PROCEDURE 

The results reported here were obtained using 
code written in C++ (with matrix multiplications 
in FORTRAN for added speed), and run on Cor- 
nell's SP-2. Gluon configurations were generated 
by the Metropolis algorithm, with 10 hits per up- 
date and 20 — 30 updates between configurations, 
which suffices to decorrelate them on the volumes 
we used. The quark propagator in each gauge 
background was found using the BiCGStab2 al- 
gorithm, supplemented by a simple type of pre- 
conditioning. 



We performed simulations at three couplings, 
/3 = Ppi — 6.8,7.1,7.4, corresponding to lattice 
spacings a — 0.40, 0.33, 0.24 fm, if one uses char- 
monium to set the scale [ ^ . We have results for 
seven quark masses at /3=6.8 and 7.1, for four at 
/3=7.4, with mp/rriT, ratios in the range between 
1.2 to 2.0. For the smaller ratios the spatial size 
of our lattices is about 2.0 fm, increasing to 2.6 — 
3.2 fm for larger ratios, except for (3— 7 A, where at 
present our largest lattices are 2.2 fm. In lattice 
units our spatial volumes are in the range 5^ to 9'^. 
The time extent of our lattices is always as least 
twice the spatial one. Since the cost of gener- 
ating gluon configurations is negligible compared 
to that of the inversion, we did not save configu- 
rations but generated independent ones for each 
quark mass. We typically produced between 200 
and 600 propagators per quark mass, the smaller 
numbers corresponding to /3=7.4. 

We calculated the time slice propagators for 
the relativistic pion, rho and nucleon operators, 
as well as for the non-relativistic nucleon, and 
the Jz = l/2 and 3/2 delta operators. Further- 
more, we obtained the pion and rho with vari- 
ous momenta. In all cases we suitably averaged 
over spin so that (in the infinite statistics limit) 
the time slice propagators for large times t are 
oc coshm(^ — t), where T is the time extension 
of the lattice. 

We used a covariant form of gaussian smear- 
ing for the quark fields at the source and sink, 

employing the operator ( 1 + e » ) ■ On 
our coarse lattices it seems that for a given /3 
one can use the same, near-optimal (in terms of 
cost/error) value of eN for all quark masses and 
hadrons. Initially we used (correlated) 2x2 ma- 
trix fits of all four combinations of sink and source 
local and smeared propagators (LL, SL, LS, SS) 
to ground and excited state cosh's. After deter- 
mining the values of eN where the SS propagator 
has only a very small excited state contribution, 
we found it more cost efficient to generate only 
the source smeared LS and SS propagators and 
fit them to a single cosh, either separately or in 
a combined fit. We used a reliable, if perhaps 
somewhat conservative method to combine effec- 
tive mass estimates from fits with different imin, 
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instead of picking the fit with the best "goodness" 
by some criterion. This is relevant in cases where 
fits from different ii^in appear to be nearly equally 
"good" . 

4. RESULTS 

We first investigated finite volume effects. We 
find that for TOp/m^ « 2, infinite volume results 
are reached within our statistical errors for the 
pion and rho at spatial sizes of around 2.0 fni, for 
the nucleon around 2.5 fm. For smaller rfiplm-^ 
ratios, somewhat smaller volumes suffice. All our 
results presented below should therefore have no 
significant finite volume errors, except perhaps 
the nucleon at A (cf. below). The delta might 
require even bigger volumes than the nucleon, and 
we will not present results for it, except to re- 
mark that the good agreement we find between 
the masses of the Jz=l/2 and 3/2 delta states, 
which are as large as am^ « 3, provide good ev- 
idence for the rotational invariance of the D234 
action on coarse lattices. 

Further such evidence is furnished by the pion 
and rho states with momenta. Defining a "speed 
of light" c via 

cV = E{pf - E{Of , (4) 

we compare in table |l| results from the D234 and 
SW actions. We see that rotational invariance, 
i.e. w 1, is restored already on the coarsest lat- 
tice for D234, in contrast to SW. This conclusion 
is even more striking for the higher momenta we 
investigated. We have also checked that c deterio- 
rates dramatically without tadpole improvement 
in the D234 action, even though it is still bet- 
ter than for SW with tadpole improvement (with 
improved glue in both cases). 

We now describe our results for the hadron 
masses. The data are shown in table ||. For nip 
we can obtain perfect linear fits to Vq + V2m^ up 
to the largest masses considered, as shown in ta- 
ble ^, and in fig. ||for /? = 7.1. Lacock and Michael 
[ proposed to use the dimensionless parameter 
J = mvdmv / dm?p at my /rap = mK* /raK = 1.8 
to characterize the vector versus pseudo-scalar 
mass relation and compare lattice studies to na- 
ture. They estimate J = 0.48(2) for the real world 



Table 1 

"Speed of light" squared (c^ in eq. |) for mesons 
with momentum p = 27r /aL, aL w 2.0 fm at 





« 0.70. 






D234 Action 


SW Action 


P 


TT p 


TT p 


6.8 


0.95(2) 0.93(3) 


0.63(2) 0.48(3) 


7.1 


0.94(3) 0.96(5) 


0.74(3) 0.55(4) 


7.4 


0.99(4) 1.00(6) 





— note that J = 0.5 if my — m?p were exactly 
constant — and J = 0.37(2) (4) as the "world av- 
erage" from quenched simulations. Our estimate 
of J from chiral fits is shown in table |^. Other es- 
timates in these proceedings are also in excellent 
agreement with ours, but do not scale on similarly 
coarse lattices. 



airip^N 




0.0 0.5 1.0 1.5 



(aTO^)2 

Figure 2. Masses of the rho (circles) and nucleon 
(diamonds) and their chiral fits at [3 = 7.1. 

We now turn to the nucleon. The masses ex- 
tracted from the relativistic and non-relativistic 
operators are found to agree well within statisti- 
cal errors. To obtain good chiral fits it is crucial, 
in general, to include a cubic term, i.e. fit mAr 
to iVo + N2ml + N^ml, cf. fig |. Our results 
are shown in table ^. For /?=7.4 we can not com- 
pletely exclude finite volume errors in the chirally 
extrapolated rriN, but they are almost certainly 
smaller than the large statistical errors we have 
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Table 2 

Light hadron masses from the D234 action 



P 


(XTTIq 




QiTftp 


QiTft ]\j 


6.8 


-0.15 


0.700(4) 


1.375(10) 


2.047(24) 




-0.08 


0.911(4) 


1.466(8) 


2.270(17) 




-0.05 


0.980(3) 


1.497(6) 


2.331(16) 







1.092(3) 


1.566(6) 


2.464(15) 




0.05 


1.206(3) 


1.638(5) 


2.598(12) 




0.1 


1.298(2) 


1.698(5) 


2.695(10) 




0.2 


1.480(3) 


1.838(5) 


2.951(14) 


7.1 


-0.14 


0.616(4) 


1.226(8) 


1.720(21) 




-0.11 


0.721(3) 


1.266(5) 


1.868(17) 




-0.06 


0.861(3) 


1.327(7) 


2.018(16) 







1.005(3) 


1.413(7) 


2.189(19) 




0.05 


1.112(2) 


1.483(5) 


2.311(15) 




0.08 


1.168(2) 


1.524(4) 


2.384(10) 




0.12 


1.251(2) 


1.581(6) 


2.474(10) 


7.4 


-0.05 


0.662(4) 


1.055(9) 


1.560(15) 







0.806(4) 


1.147(7) 


1.739(10) 




0.05 


0.921(5) 


1.224(9) 


1.869(14) 




0.1 


1.031(3) 


1.304(4) 


2.010(9) 



Table 3 

Chiral fit parameters (in lattice units) of the rho 



mass 







V2 


iVDF 


Q 


6.8 


1.238(7) 


0.274(5) 


5 


0.98 


7.1 


1.108(6) 


0.303(5) 


5 


0.94 


7.4 


0.889(13) 


0.390(14) 


2 


0.67 



in this case. In table ^ we show our estimates of 
mN I'Tfip (=1.22 in nature) in the chiral limit. 

So far the results described for the D234 ac- 
tion have been spectacular. There is however one 
problem, which is revealed in table |[ the scale 
determined from the chirally extrapolated p mass 
(set to 770 MeV) does not agree with that from 
the charmonium IS— IP splitting (00), which is 
known to scale quite well [ ^ . In view of the ex- 
cellent properties of (? and J we interpret this 
as indicating a small additive shift of all masses 
(vanishing, of course, as a ^ 0) that depends, it 
seems, only rather weakly on the momentum or 
type of the particle, at least for mesons. We are 
investigating whether the doublers are responsi- 



Table 4 

Chiral fit parameters (in lattice units) of the nu- 
cleon mass 





N2 


A^s iVDF Q 


6.8 1.66(8) 
7.1 1.35(5) 
7.4 1.09(11) 


0.99(21) 
1.30(18) 
1.44(46) 


-0.29(13) 3 
-0.47(12) 4 
-0.56(35) 1 


0.91 
0.80 
0.55 


Table 5 

Selected quantities relevant for scaling, mj^jmp 
is extrapolated to zero pion mass, is from 


(5 J 


niM/mp 


(MeV) 


a-ip/ap 


6.8 0.384(8) 
7.1 0.381(8) 
7.4 0.395(18) 


1.34(6) 
1.22(5) 
1.23(12) 


497(4) 
606(5) 
821(10) 


1.251(12) 
1.147(11) 
1.055(20) 



ble for this shift. There are various strategies 
to "push away" the doublers, and we are hopeful 
that a completely satisfactory fermionic action for 
coarse lattices can be found. 
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